Electron correlation induced transverse delocalization and longitudinal confinement in 
excited states of phenyl-substituted polyacetylenes 
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Electron-electron interactions in general lead to both ground state and excited state confinement. 
We show, however, that in phenyl-substituted polyacetylenes electron-electron interactions cause 
enhanced delocalization of quasiparticles in the optically excited state horn the backbone polyene 
chain into the phenyl groups, which in turn leads to enhanced confinement in the chain direction. 
This co-operative delocalization-confinement lowers the energy of the one-photon state and raises 
the relative energy of the lowest two-photon state. The two-photon state is slightly below the optical 
state in mono-phenyl substituted polyacetylenes, but above the optical state in di-phenyl substituted 
polyacetylenes, thereby explaining the strong photoluminescence of the latter class of materials. We 
present a detailed mechanism of the crossover in the energies of the one- and two-photon states 
in these systems. In addition, we calculate the optical absorption spectra over a wide wavelength 
region, and make specific predictions for the polarizations of low and high energy transitions that 
can be tested on oriented samples. Within existing theories of light emission from 7r-conjugated 
polymers, strong photoluminescence should be restricted to materials whose optical gaps are larger 
than that of trans-polyacetylene. The present work show that conceptually at least, it is possible to 
have light emission from systems with smaller optical gaps. 

42.70.Jk,71.20.Rv,71.35.-y,78.30.Jw 

I. INTRODUCTION 



The absence of photoluminescence (PL) in linear 
polyenes, trans-polyacetylene (t-PA) and the polydi- 
acetylenes (PDAs) is a well understood electron correla- 
tion effect. Excited states in these centrosymmetric sys- 
tems have Kg or B„ symmetry, which can be further clas- 
sified according to their charge- conjugation symmetryEJ. 
The dipole selection rule allows optical transitions only 
between Ag and B„ states of opposite charge-conjugation 
symmetries, and the optical transition from the lAg 
ground state to the 1B„ state with opposite charge- 
conjugation symmetry is strongly allowed. Due to the 
moderate strength of the on-site electron-electron (e-e) 
interaction (the Hubbard U) in these systems, hcMvgver, 




the lowest dipole-forbidden 2kg state occurs belowErH the 
1B„, and the optically pumped 1B„ decays in ultrafast 
times to the 2Ag, radiative transition from which to the 
lAg is forbidden. 

Strong PL from 7r-conjugated polymers like poly(para- 
phenylene) (PPP), poly(para-phenylenevinylene) (PPV) 
and their derivatives is also by now„syell understood. 
As shown by Soos and coUaboratordaa the lowest ex- 
cited states of these systems can be understood within 
effective linear chain models with large bond alterna- 
tion. Schematically, materials like PPP and PPV can be 
thought of as obtained by replacement of all or alternate 
double bonds in t-PA with benzene moieties, a process 
we term as "bond substitution", that in effect increases 
the bond alternation of the effective linear chainQ. The 
enhanced bond alternation leads to a reversed energy or- 
dering E(2Ag) > E(1B„) (where E(...) is the energy of 
the state in question), which is conducive to strong light 
emission. 



(a) (b) 

FIG. 1. The unit cells of (a) PPA and (b) PDPA. The 
phenyl rings are rotated with respect to the y-axis, which is 
transverse to the axis of the polyene backbone (a;-axis) 

The recent observation of strong PL in phen¥j-|-dis- 
ubstituted polyacetylenes (PDPAs) (see Fig. pjlB is 
difficult to explain within the above simple theoretical 
picture. These systems consist of a backbone polyene 
chain, the hydrogen-atoms of which have been replaced 
with phenyl groupalj, a process that we term as "site 
substitution" , in order to distinguish it from bond sub- 
stitution. Since the Hubbard U for the backbone polyene 
carbon atoms (the repulsion between two 7r-electrons oc- 
cupying the same backbone atom), is an atomic prop- 
erty, it should have the same magnitude in the PDPAs 
and polyenes. Since the one-electron hopping integrals 
between the backbone atoms are also expected to have 
similar magnitudes, the observed PL is perplexing. Sev- 
eral different explanations of this apparently peculiar be- 
havior are possible, and one goal of the present paper is 
to arrive at the correct mechanism of the PL. A bsief pre- 
sentation of our work has been reported earlierlla. Here 
we present more detailed work with a broader goal, viz., 
to understand the 2Ag - 1B„ crossover that we show oc- 
curs in these systems from a global perspective for the 
family of vr-conjugated polymers. Specifically, we are in- 
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terested in the mechanism of the crossover, and not just 
the result. We show that there are two different man- 
ifestations of the so-called "electron correlation effect", 
and as far as the relative location of the 2Ag is con- 
cerned, the two effects are competing. The effect of the 
spin-dependent Hubbard correlation is to always lower 
the energy of the 2Ag and to increase the energy of the 
IBu. A second correlation effect can be the confinement 
of electrons and holes in the excited IB^ state, rfp^h has 
been discussed only in the context of excitonstSllZl. Al- 
though the Hubbard U causes ground state localization, 
it has a weak effect on the electron-hole confinement in 
the excited state, and the origin of any such confinement 
therefore has to be different. There can be multiple ori- 
gins of this confinement, but in perhaps all these cases 
the consequence is either the lowering of E(lB,j), or the 
raising of E(2Ag), or both, i.e., the opposite of the effect 
due to the Hubbard U . Our aim in the present paper 
is to give a clear understanding of how this confinement 
occurs in the PDPAs. 

A second goal of the present work is to understand 
and predict the optical absorption in phenyl-substituted 
acetylenes over a wide wavelength region. In the past, 
theoretical and experimental work on the polarizations of 
the high energy absorptions in oriented samples of PPV 
have considerably increased our und|eijSjta|nding of elec- 
tron correlation effects in this systcmllrEil, and it is our 
hope that similar polarized absorption measurements will 
be carried out on oriented samples of phenyl acetylenes 
as they become available. 

In order to discuss the applications of our results to 
real materials, we need to clarify several aspects of the 
experiments with the phenyl-substituted polyacetylenes. 
Before proceeding further, we therefore-.aresgnt a brief 
overview of the experimental situationlalj'E3E3. Three 
different classes of phenyl acetylenes have been synthe- 
sized, which include, (i) poly-phenylacetylenes (PPAs), 
which are obtained by substituting alternate hydro- 
gen atoms of t-PA with phenyl groupsllj, (ii) poly- 
phenylalkylpolyacetylenes (PAPAs), in which alternate 
hydrogens of t-PA are replaced with alkyl and phenyl 
groupsM, respectively, and (iii) the PDPAs, already men- 
tioned above. In Figs. |l|(a) and (b) we have shown the 
structures of PPA and PDPA, respectively. The PPAs 
are weak emitters, while whether or not the PAPAs emit 
light depends on the size of the alkyl group, - for small 
alkyl groups the PL quantum efficiency is small while 
for large alkyl groups this can be as large as in PDPAs 
(see Fig. 6 in reference In contrast, all PDPAs are 
strong emitters of light, in the green or blue region of the 
spectrum. Raman scattering studies that have probed 
the backbone C = C bonds have determined that— the 
conjugation lengths of PAPAs and PDPAs are shortEl (5 
- 7 double bonds along the backbone according to the 
authors). The finite conjugation length is perhaps to be 
anticipated, based on the large sidegroups that are apt 
to lead to twists and bends involving the backbone main 
chain. Finally, in both PPAs and PDPAs (see Fig. |^) 



steric effects involving the phenyl groups would be con- 
siderable, and moderate phenyl group rotations are to be 
expected. There is currently no experimental informa- 
tion on the extent of this ring torsion. 

In principle, there can be multiple reasons for the ob- 
served light emission in PDPAs. First, the light emission 
can be a simple consequence of the loss of C2h symmetry 
due to ring torsions. Within this scenario the nominally 
2Ag occurs below the 1B„, but in the absence of strict 
symmetry classification, this state is emissivellJ. Second, 
the energy ordering is reversed, i.e., E(2Ag) > E(1B„), 
but this is a simple consequence of finite conjugation 
lengths, - it is, for example, believed that in the short- 
est polyenes the 2Ag might occur above the 1B„. Third, 
it is conceivable that the light emission is a consequence 
of an intricate electronic mechanism that involves the 
electron hopping among the phenyl groups. Referring to 
Fig. 1^, for example, even with ring torsions the distances 
between consecutive phenyl groups on the same side of 
the backbone polyene chain in both PPA and PDPA are 
rather small. If these are now nearly coplanar (i.e., if 
the ring torsions are ordered) there should be consider- 
able overlaps between the 7r-electron clouds of the piienyl 
groups, leading to electron hoppings between themE3, as 
is well-known in the literatiire on segregated stack or- 
ganic charge-transfer solidsa. Fourth, it is conceivable 
that phenyl-substitution leads to enhanced bond alter- 
nation on the backbone polyene chain, and this is what 
leads to the reversed energy_prdering and light emission, 
exactly as in PPV and PPPo. Finally, light emission can 
be a consequence of confinement or localization along the 
backbone, which is accompanied, by enhanced delocaliza- 
tion into the phenyl sidegroupslij. We have ignored in the 
above the possibility that excited state ordering E(2Ag) 
> E(IB„) is due to the relaxation of the 1B„ exciton, 
as the 1B„ - 2Ag energy gap in t-PA is known to be 
fairly large, and even if exciton relaxation plays a role, 
this would still require that this energy gap in PDPAs 
has shrunk due to an electronic mechanism. 

We point out that the experimental results already pre- 
clude most of the possible scenarios. The loss of G2h 
symmetry is stronger in PPAs than in PDPAs, and if 
the light emission was from the nominally 2Ag state, one 
would have expected these to be more strongly emissive 
than the PDPAs. Furthermore, in our explicit calcula- 
tions for PPA, we have found that the transition dipole 
coupling between the nominally 2Ag and the nominally 
lAg is nonzero, but its magnitude remains very small (we 
therefore continue to refer to eigenstates as Ag and B„, 
even though this classification is not strictly valid). Lack 
of C2h symmetry is therefore not the origin of the strong 
PL. Finite conjugation length alone cannot explain, the 
light emission, as the estimated conjugation lengthsEy are 
in the regiop-|Where the 2Ag occurs below the 1B„ in 
the polyenesE3. The direct electron hopping between the 
phenyl groups do introduce a novel conjugation channel 
absent in other 7r-conjugated polymers, but this effect 
should be comparable in the PPAs, PAPAs and PDPAs. 
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The absence of strong PL in the PPAs and many PAPAs 
then again precludes this mechanism. Indeed, as we show 
from exphcit calculations, this particular interaction af- 
fects the lowest energy allowed and forbidden optical gaps 
very weakly, although it can have a strong effect on the 
optical absorption in the high energy region. Thus from 
the experiments alone, we believe that there arc only 
two possible explanations of the strong PL in PDPAs, 
viz., enhanced bond alternation and an electron correla- 
tion induced confinement-delocalization effect. Our pre- 
vious work had already shown that the bond alternation 
is smaller in the PDPAs than in t-PA. In the present work 
we show further detailed work, emphasizing on the mech- 
anism behind the reduced bond alternation, rather than 
the result itself. Our aim is to show that this effect and 
the confinement effect are ultimately related, and even 
though the 2Ag - 1B„ crossover is a many-body prob- 
lem, the final results could have been perhaps anticipated 
from a combination of our understanding of one-electron 
theory and confinement effects due to e-e interactions. 

In the rest of the paper we present detailed calcula- 
tions for PPA and PDPA oligomers and polymers. We 
do not present explicit calculations for PAPA, the un- 
derstanding being that the effect of the alkyl groups on 
the electronic structure is by and large the same as that 
of the hydrogen atoms, with the only difference that the 
site energies (electronegativities) of the backbone carbon 
atoms to which the alkyl groups are bonded might be 
weakly different. We postpone comments on the strong 
PL from PAPAs with very large alkyl groups until the 
last section of the paper. In section II we present our 
Hamiltonian, along with a detailed discussion of the pa- 
rameters. In section III we give a brief discussion of 
known results concerning the effect of confinement on the 
2Ag - 1B„ crossover in linear chains. These discussions 
point out already the necessity to have a precise under- 
standing of the electronic structures of PPA and PDPA 
within the one-electron limit of our Hamiltonian, our re- 
sults for which are presented in section IV. In section V, 
we present details of the optical absorption within the 
correlated Hamiltonian within the singles-configuration 
interaction (SCI) approach. The goal of this section is 
two-fold. First, we show that the electron-hole confine- 
ment along the backbone has a synergistic relationship 
with the delocalization into the phenyl rings. A second 
goal is to make specific predictions that can be tested 
out in absorption studies of oriented samples in the fu- 
ture. In section VI, we present the results of multiple 
reference singles afid doubles configuration interaction 
(MRSDCI) studies^ on short PPA and PDPA oligomers 
to prove the reversed energy ordering in these systems. 
These calculations are considerably morCpSophisticated 
than our previously reported calculationscJ, but never- 
theless, slightly longer oligomers could be investigated 
here. In section VII we present our concluding discus- 
sions, emphasizing the fact that the mechanism of 2Ag - 
1B„ crossover found here introduces the novel possibil- 
ity of synthesizing light emitting 7r-conjugated polymers 



with small optical gap; 
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II. THEORETICAL MODEL 

Our calculations are within the 7r-electron model for 
PPAs and PDPAs, as the ordering of the lowest energy 
states depends primarily on the electron correlation ef- 
fects among the 7r-electrons. Even with this restriction, 
the complete theoretical model for PPAs and PDPAs is 
rather complicated, as this must include the effects of 
the electron-phonon (e-ph) couplings along the backbone 
polyene chain, electron-electron (e-e) interactions among 
all the TT-electrons, and all effects of phenyl-group rota- 
tions. Phenyl-group rotations modify the one-electron 
hoppings between the phenyl groups and the backbone 
chain and the intersite Coulomb interactions involving 
the carbon atoms of the phenyl groups, and in addition, 
introduce cr-typc bondings between the 7r-molecular or- 
bitals (MOs) of neighboring phenyl groups on the same 
side of the polyene chain. This last interaction is maxi- 
mized when the rotational angle is ordered and 90°, but 
such large rotation would completely destroy the phenyl- 
polyene conjugation, leading to an effectively isolated 
polyene chain and isolated benzene molecules. Since the 
properties of PPAs and PDPAs are different from those of 
isolated polyene chains, the rotation angle must be con- 
siderably less than 90°, in which case the one-electron 
hoppings between neighboring phenyl groups can involve 
only a few of the carbon atoms of each phenyl group (see 
Fig. 0). 

From the above, the overall Hamiltonian for PPAs and 
PDPAs is written as. 



H — Hn + Hfi + Hi 
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Here He and Hb are the one-electron Hamiltonians for 
the backbone chain carbons and the benzene units, re- 
spectively, HcB and Mbb are the one-electron hoppings 
between the chain and the phenyl units, and between the 
phenyl groups themselves, respectively, and Hee is the e-e 
interaction. The individual terms can now be written as. 
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In the above, k, k' are carbon atoms on the polyene 
backbone, /i, v are carbon atoms belonging to the phenyl 
groups, M is a composite site consisting of a phenyl group 
and a polyene carbon, (...) implies nearest neighbors, 

and B.j.MM' = Y.ai4.M,a^3-M' ■'y + ^-c-)- The io, i± 
and tcr(/i, J^) are matrix elements corresponding to one- 
electron hops, where t„ depends on the particular /i, v 
being considered. In Hq, a is the e-ph coupling constant, 
K the spring constant, and /^k,M = (ufc+i,j\/+i - Uk,M), 
where Uk,Ai is the displacement of the Mth composite site 
from equilibrium. In Hcb, the sum over fi is restricted 
to atoms of the phenyl groups that are directly bonded 
to backbone carbon atoms. The prime on the summa- 
tion in Hbb indicates that the sum is restricted to a few 
carbon atoms of neighboring phenyl groups that are on 
the same side of the polyene chain (i.e., M' = M ± 2), 
the actual number depending upon the extent of phenyl 
group rotation. H^,, is the e-c interaction, with i and M 
including now all atoms. The Coulomb interactions are 
parametrized according to the Ohno relationshipc3, 

V^^J,M,N = C//(l + Q.&inRl,^MMY'^ , (2) 

where U = 11.13 eV and Rij,M.N is the distance in A be- 
tween the zth carbon in unit M and the jth carbon in 
unit N . The Ohno parametrization is for small molecules, 
and it is conceivable that the Coulomb parameters for 
the polymeric samples are somewhat smallerE3. Our aim, 
however, is to obtain a qualitative understanding of the 
role played by e-e interactions on the optical properties 
of the system. We have therefore employed large H^e 
to show that in spite of these large interactions energy 
ordering E(2Ag) > E(1B„) will occur in the substituted 
polyacetylenes. The same considerations apply to the 
other parameters (see below), viz., the actual numerical 
values are relatively unimportant, as our goal is to ex- 
tract the essential physics of these systems. 

In our calculations the backbone polyene atoms were 
assumed to lie in the xy plane, with the x direction de- 
fined as the polymer axis (see Fig. 1). Due to steric 
effects, the phenyl rings will not be in the xy plane, but 
will be rotated with respect to the y axis. All intra- 
phenyl bond lengths are taken to be 1.4 A and the bond 
between the backbone atoms and the phenyl groups is 
taken to be a true single bond with length 1.54 A. In our 
calculations of bond-alternation along the polyene back- 
bone (section IV), the backbone polyene bond lengths 
are assumed to be 1.4 A in the undistorted configuration, 
while the bond lengths are calculated self-consistently in 
the distorted configuration. The calculations of optical 
absorption are for finite oligomers, in which we chose the 
backbone bond lengths the same as in t-PA, viz., 1.45 A 



and 1.35 A, for reasons to be discussed later. We con- 
sidered both parallel and antiparallel ordered rotations of 
the phenyl groups bonded to the nearest neighbor carbon 
atoms in PDPA (note that these occur on opposite sides 
of the backbone, the phenyl groups on the carbon atoms 
on the same side of the chain are parallel in both cases) 
with rotation angles of 30° (this value of the torsional 
angle is of course arbitrary; on the other hand torsional 
angles larger than 45° will decrease t± drastically). How- 
ever, the differences in numerical results (energies as well 
as transition dipole couplings) for the two cases were in- 
significant. Therefore, we restricted our studies to the 
parallel configuration for the phenyl groups, rendering 
the symmetry group of PDPA to be Ci, with identity 
and inversion operators as symmetry elements. PPA, on 
the other hand, does not possess any point group sym- 
metry. 

As far as the hopping matrix elements are concerned, 
we assumed to = 2.4 eV. Along the backbone carbon 
atoms, whenever rigid bond alternation was assumed 
along the chain, standard values of 2.2 eV and 2.6 eV 
were chosen for the long and the short bonds, respec- 
tively. The aforesaid rotation of the phenyl groups with 
respect to the y axis will have the effect of reducing the 
hopping integral tj_ connecting the phenyl rings to the 
t-PA backbone, as compared to the standard value of 2.4 
eV. Various possible values for t± were considered in our 
calculations. We postpone the discussion of our choice of 
to-(/i, ly) until later. 

III. CONFINEMENT AND THE 2A9-1B„ 
CROSSOVER 

Although the theory of 2Ag - IBu crossover is well 
understood in linear chains, it is useful to present a brief 
review in view of what follows. Consider the rigid bond 
limit of Eq. (1) for the case of linear chains, with t± = 
and the chain hopping integrals to{l±S). In the limit of 
large U, the ground state has all atoms singly occupied by 
electrons and the 2Ag is a spin excitation that is sum of 
two triplets, each localized on the double bonds, and has 
the energy 4i^(l -I- / [U — V), where V is the Coulomb 
interaction between nearest neighbor atoms on the linear 
chainEra. 

For finite U , the overall wavefunction is a superposi- 
tion of this spin excitation and a charge excitation, and 
the 2Ag occurs below the 1B„ if the contribution of the 
spin excitation configuration dominates over that of the 
charge excitation configuration. From configuration in- 
teraction principles, the relative contribution of the spin 
excitation decreases as its energy increases, and therefore 
both increasing 5 (for fixed U and V) and increasing V 
(for fixed U and S) will increase the energy of the 2Ag 
relative to the 1B„. This can be understood physically in 
the limit of large 6 orV rather easily. The spin excitation 
that is relevant in the description of the 2Ag is obtained 
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FIG. 2. Band structures of (a) PPA and (b) PDPA for t± = lA eV and U = 0, and (c) PPA for tj_=lA eV and = 1.0 
eV. Here thick solid lines represent polyene chain-derived bands, with zero bond alternation in (a) and (b), and with finite 
bond alternation in (c). Finite bond alternation in (a) and (b) opens a gap at the Fermi level, and the highest valence and 
lowest conduction bands in these cases are shown as thin dashed-dotted lines. The thin solid lines represent benzene-derived 
localized bands (see text) and the thick dashed-dotted lines in (b) represent benzene-derived delocalized bands. There occur 
delocalized benzene-derived bands in (a) which are nearly degenerate with the flat localized bands and are not visible within 
our resolution. 



from the localized ground state by two consecutive elec- 
tron hops in opposite directions, the first creating a vir- 
tual state with a doubly occupied site and a vacancy, 
and the second destroying this. Now increasing 6 or V 
both have the tendency to confine the double occupancy 
and the vacancy in the virtual state to nearest neighbor 
dimer cthylenic units, and since the dimer does not have 
a spin singlet covalent excitation, there always exist a 6 
or V where E(2Ag) > E(1B„), with the only condition 
that S is nonzero to begin with. Thus any interaction 
that confines the double occupancy and the vacancy in 
the virtual state should increase the relative energy of 
the 2Ag. We refer to references |[|| for further details on 
the effect of S, but point out a fundamental difference 
between the effect of S and V on the IB^,. While increas- 
ing d (for fixed U and V) increases E(1B„), increasing V 
(for fixed U and S) decreases E(1B„) (since the correla- 
tion contribution to the optical gap is linear in U — V). 
Although increasing V is not a practical proposition in 
linear chains, the more important point is that in princi- 
ple at least, energy ordering conducive to light emission 
does not necessarily require that E(1B„) be larger than 
in the polyenes or t-PA. 

Recognition that any interaction that promotes con- 
finement of the vacancy and the double occupancy in the 
virtual state allows us to understand the mechanism of 
a possible 2Ag - 1B„ crossover in PDPA by probing the 
1B„ in detail, since the same interaction will also lead 
to confinement in the 1B„. Now, probing the 1B„ the- 
oretically is much simpler than probing the 2Ag. This 
is because the SCI approximation, which is easy to im- 
plement, is known to give a reasonable description of the 
1B„, but not that of the 2Ag. Any possible confinement 
in the 1B„ in PPA and PDPA can only occur due to 
conjugation with the phenyl sidegroups, and in order to 
understand the effects of this conjugation in detail we 



have to start from one-electron theory, as is done in the 
next section. 



IV. ONE-ELECTRON THEORY 

In this section we present the results of our calcula- 
tions of band structures, bond alternations, and optical 
absorptions in PPA and PDPA based upon one-electron 
theory [Hee — within Eq. (1)). It will be seen that 
the difference between PPA and PDPA on the one hand, 
and polyenes on the other, can be largely understood al- 
ready within one-electron theory, and the effects of e-e 
interactions can be anticipated based upon these results. 



A. Band Structure 

The band structures were calculated within the rigid 
bond approximation, and two different cases, viz., equal 
bond lengths and bond lengths corresponding to those of 
t-PA along the polyene backbone were considered. Fur- 
thermore, since the number of carbon atoms to be re- 
tained in Hbb is somewhat arbitrary (as the exact ex- 
tent of the overlap between the tt-MOs of neighboring 
phenyl groups is difficult to estimate) we perform two 
different sets of band structure calculations, with Hbb 
both zero and nonzero. As these calculations indicate, 
the strongest effects of phenyl substitution in the region 
close to the Fermi level are those due to Hcb, with Hbb 
playing a smaller quantitative role. Hbb, however, plays 
a stronger role away from the Fermi level. 

For Hbb — the one-electron part of Hamiltonian 
(1) is the simple Hiickel model. In this case, the con- 
sequence of nonzero Hcb is simple hybridization be- 
tween the polyene bands and benzene molecular orbitals 
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(MOs). The highest (lowest) occupied (unoccupied) ben- 
zene MOs (HOMOs and LUMOs) are doubly degenerate, 
with one member of each set of doubly degenerate lev- 
els being delocalized over all the carbon atoms of the 
benzene molecule and the other localized over four of 
the six carbonscS. Henceforth, we refer to the delocal- 
ized HOMO (LUMO) as d (d*) and the locahzed HOMO 
(LUMO) as / (/*). We expect very strong hybridization 
between the polyene valence (conduction) bands and the 
d (d*) MOs, but the I and /* MOs are-estp^ed to retain 
their characters, as in the case of PPVE3Ejc3'E2l. In addi- 
tion to the d, d* frontier MOs, each benzene molecule also 
has a delocalized low (high) energy occupied (unoccu- 
pied) MO, hybridization between which and the polyene- 
derived bands is expected to be smaller (except at the 
edge of the Brillouin zone) due to the large one-electron 
energy difference. 

In Figs. I (a) and (b) we have shown the band struc- 
tures for PPA and PDPA, respectively, for Hbb — 0, 
but nonzero Hcb, with t± = 1.4 eV. We have verified 
that the band structure remains qualitatively the same 
for all t± between 1 eV and 2.0 eV. For the case of zero 
bond alternation, there exists a Fermi surface degeneracy, 
exactly as in the case of t-PA, indicating unconditional 
bond alternation in the infinite chain limit. The band 
structures are nearly identical for zero and nonzero bond 
alternation, except at the center of the Brillouin zone, 
where a gap opens up for nonzero bond alternation in 
both PPA and PDPA. 

For phenyl group rotations of 30° - 45°, Hbb has sig- 
nificant contributions only from the overlaps between the 
carbon atoms 2 and 6', and between atoms 3 and 5' in 
Fig. 1^ (here 6' and 5' refer to carbon atoms equivalent to 
6 and 5, on the next phenyl group on the same side of the 
polyene chain) . We will therefore ignore the more distant 
to-(/i, ly) and will henceforth refer to the nearest neighbor 
terms that we retain as simply t^. In the absence of 
theoretical and experimental information concerning the 
magnitudes of these hopping integrals, we have chosen 
relatively large values of 1.0 eV and 2.0 eV, to determine 
the largest effects of Hbb- The band structure of PPA 
for the case of = 1.0 eV is shown in Fig. ||(c). There 
are two effects of nonzero Hb b ■ First, charge conjugation 
symmetry is broken. Second, there is now much larger 
dispersion in the "benzene-derived" bands. However, the 
gap between the highest valence band and the lowest con- 
duction band remains largely unaffected, a result that 
remains true even for = 2.0 eV (not shown). Due 
to the doubling of the unit cell upon inclusion of Hbb 
the Brillouin zone is halved and the number of bands for 
PPA is twice those in Fig. |(a). In the case of PDPA 
this number is larger by almost another factor of two, 
and we have therefore not shown the band structure for 
PDPA for nonzero ta. As in PPA, however, Hbb has 
an insignificant effect on the energy gaps in the low en- 
ergy region. We therefore do not expect the low energy 
physics (especially the energy ordering upon inclusion of 
Hee) to be altered by Hbb- On the other hand, because 



of the larger dispersion in the benzene-derived bands, 
nonzero Hbb leads to chain-to-benzene charge-transfer 
absorptions at relatively lower energy (see below). 

From the band structures alone, the strong hybridiza- 
tion that occurs between the polyene bands and the ben- 
zene delocalized HOMO and LUMO are obvious only 
away from the center of the reduced Brillouin zone (k — 
7r/2a in Figs. 2(a) and (b) and fc = in Fig. 2(c)). A 
more quantitative measure of the hybridization, however, 
is the contribution of the chain atoms to the band wave- 
functions. We have calculated this quantity for the one- 
electron wavefunction corresponding to the Fermi level 
kp, defined as, p^k^ = J2techain Cfc^-.iCfcf where Cfc^,^ is 
the coefficient of the ith chain atom to the wavefunction 
a.t k = kp, for both PPA and PDPA as a function of t±. 
In Fig. ^a) we have summarized our results, where we 
see that Pcfep decreases rapidly with t±. Similar calcu- 
lations were done also for one-electron levels away from 
kp. In general, the chain contribution to the highest oc- 
cupied and lowest unoccupied bands decreases with in- 
creasing absolute value of the one-electron energy (there 
of course exists a sum rule, which implies that the nom- 
inally benzene-derived bands away from the Fermi level 
have considerable contribution from the polyene chain 
atoms). Since the Peierls bond-alternation in t-PA is a 
strictly one-dimensional effect, the large decrease in Pc.kp 
with tj_ indicates that the bond-alternation in PPA and 
PDPA should be smaller than in t-PA. 



B. Bond Alternation 

For bond alternation within a periodic configuration, 
Uk,M — (~l)^^uo in Eq. (1). We have calculated the op- 
timal bond alternation uq, within the one-electron limit 
{Hee = 0) of Eq. (1). The magnitude of uq depends on 
the dimensionless e-ph coupling constant a^/Kto, but 
the qualitative behavior as a function of t± does not. We. 
have therefore chosen the original Su-Schrieffer-Heegert3 
values of a (4.1 eV/A) and K (21 eV/A^), for direct 
comparisons to previous one-electron theories of t-PA. In 
reference |l^, we had shown the behavior of finite PDPA 
rings with variable size for a few fixed t±. Here we show 
the behavior of both PPA and PDPA for < < 2.0 
eV in Fig. ||(b) . These results are for periodic rings con- 
taining 100 backbone carbon atoms, where the optimal 
Uq has reached its saturation value in all cases. The close 
relationship between the optimal uq and p^kp is obvious 
from the natures of these curves. 

As has been argued in reference |l^, nonzero H^e will 
change the magnitude of uq in each polymer, but will 
not change the trend seen in Fig. ^ (b), viz., a reduction 
in the bond alternation with increasing t± . We therefore 
conclude that the observed PL in PPA and PDPA cannot 
be a consequence of large effective bond alternation, as is 
true in the bond-substituted materials PPP and PPVs. 
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FIG. 3. (a) The contributions by the polyene chain atoms to the one-electron wavefunction ai k = kp in PPA and PDPA, as 
a function of tj_ for ta = 0. Nonzero to- has almost no effect on pc,kp- (b) The self-consistent bond alternation (in nanometers) 
within the one-electron limit of Eq. (1), as a function of t±. The continuous curves are for ta = 0, while the stars (★) and 
circles (o) correspond to fo = 1.0 eV for PDPA and PPA, respectively, again indicating the insignificant effect of ta on the low 
energy properties. 



C. Optical absorption 

Although our interest lies primarily in calculating the 
optical absorption with nonzero i?ee, it is useful to dis- 
cuss the i?ee — limit first. In the case of PPV, compar- 
isons of absorption spectra calculated within one-electron 
and many-electron theories on the jona. hand, and with 
experimental spectra on the otherOo, have led to a 
clear understanding of the effects of the e-e interaction. 
Although oriented samples of phenyl-substituted poly- 
acetylenes do not exist at the moment, we give details of 
the calculated polarizations here, in anticipation of future 
experimental work. 

The optical absorption within one-electron theory can 
be anticipated from the structures of PPA and PDPA. In 
the absence of Re b and Hbb we expect isolated low en- 
ergy polyene chain absorption polarized primarily along 
the X-axis, and degenerate benzene transitions at high en- 
ergy which can be characterized as d — > d* (j/-polarized) , 
I — > I* (y-polarized) and d ^ I* and I d* {x- 
polarized) (these polarizations can be anticipated from 
previous work on PPV, as the unit cell in PDPA is sim- 
ply the trans-stilbene molecule). The anticipated effects 
of nonzero Hqb are: (i) enhanced j/-character in the low- 
est energy chain absorption due to hybridization with the 
benzene delocalized MOs, (ii) splitting of the high energy 
molecular benzene transitions, and (iii) new degenerate 
chain — > benzene and benzene — > chain charge-transfer 
type absorptions at intermediate energies. The charge- 
transfer type absorptions will include transitions from 
and to benzene d, d* and I* orbitals. Of these, the for- 
mer (involving d, d* MOs) are expected to have mixed 
polarization (since these are hybridized with the chain- 
derived bands that have some y-character, see above), 
but the latter (involving I, I* MOs) should be primar- 
ily x-polarized. From the band structures in Figs. 0(c), 



for nonzero Hbb "we expect even larger splitting in the 
molecular transitions, and thus the width of the absorp- 
tion band in this region should increase. Due to broken 
charge-conjugation symmetry, the chain -—>■ benzene and 
benzene — > chain absorptions are also nondegenerate, and 
therefore the absorption band in the intermediate energy 
should also be broader. 

In Figs. U (a) and (b) we have shown the absorp- 
tion spectra of finite PPA and PDPA oligomers with 
eight double bonds along the polyene chain for nonzero 
HcB {tj_ = 1.4 eV) but Hbb = 0. Our choice of fi- 
nite oligomers is based on, the finite conjugation lengths 
of the existing materialalj. The hopping integrals along 
the polyene backbone here are 2.2 eV and 2.6 eV, cor- 
responding to those of t-PA. This is primarily because 
it is difficult to estimate the extent of the bond alterna- 
tions in finite chains, but also because the polarizations 
and the oscillator strengths of the various transitions de- 
pend weakly on the magnitudes of the hopping integrals, 
and even the energies of the higher energy transitions do 
not depend on these, as we have verified from explicit 
calculations. We have indicated in Fig. 4 the polariza- 
tions of the different absorption bands. The oscillator 
strength of the lowest energy chain — > chain absorption 
is the same in both cases, but that of the high energy 
molecular absorption (centered around 2tQ — 4.8 eV) in 
PDPA is nearly twice that in PPA. The large width of this 
absorption is due to the nondegeneracy of the molecular 
transitions. At intermediate energy, the weak second and 
third peaks correspond to chain HOMO - 1 ^ LUMO -I- 
1 and HOMO - 2 ^ LUMO + 2 transitions, following 
which there occur the charge-transfer transitions. The 
absorption spectra for other t± can be guessed from Fig. 
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FIG. 4. The calculated x (solid line) and y (dotted lines) components of the optical absorption spectra within one-electron 
theory for finite PPA (Figs, (a) and (c)) and PDPA (Figs, (b) and (d)) oligomers with eight double bonds on the backbone 
polyene chain. The y-axis scales are different for PPA and PDPA. In all cases t± = 1.4 eV, while to- = in Figs, (a) and (b) 
and 1.0 eV in Figs, (c) and (d)). 



^(a) and (b). Specifically, larger (smaller) the t_L, the 
broader (narrower) is the molecular absorption and the 
stronger (weaker) are the charge-transfer type absorp- 
tions. 

The absorption spectra for the case of nonzero Hbb 
are shown in Figs, ^c) and (d). With increasing the 
molecular absorptions are broader, the charge-transfer 
transitions are redshifted and acquire greater width. 



V. OPTICAL ABSORPTION FOR NONZERO 



sonable descriptions of the odd parity B„ states (espe- 
cially the lowest i?„ state, the 1B„), even though it is a 
poor approximation for the even parity Kg states. In Ta- 
ble I we have given the optical gaps of PPA and PDPA 
oligomers with eight double bonds for Hcb 0, but 
ilsB = (reference |l^ discusses the case of five doubel 
bonds) . These calculations are within the rigid bond ap- 
proximation, with the hopping integrals along the chain 
2.2 eV and 2.6 eV, for single and double bonds, respec- 
tively. We show later that as within one-electron theory, 
the role played by Hbb in the low energy region is weak. 



There are two goals of our work in this section. First, 
we demonstrate that the polarization characters of the 
lowest absorption in PPA and PDPA change radically 
for nonzero iJee, and that this change in the polarization 
character provides the hint to understanding the mech- 
anism of the 2Ag - 1B„ crossover that occurs in PDPA. 
Second, a strong x-polarized charge-transfer type absorp- 
tion appears at intermediate energy when both H^e and 
to- are nonzero. This intermediate energy absorption may 
have already been seen experimentally in unoriented sam- 
ples (see Fig. 2 in reference ^. In the following we discuss 
the energies and the polarizations of the lowest energy 
absorption and the overall absorption spectra separately. 



TABLE I. Optical gaps of oligomers of PDPA and PPA, 
containing eight double bonds in the backbone chain, within 
the Hiickel and the Coulomb correlated models, as a function 
of t±. Calculations based upon the correlated model were 
performed within the SCI approximation. Note that the case 
of fx ~ 0.0 corresponds to the eight double-bond polyene. All 
energies are in the units of eV. 

tx PPA PDPA 





Hiickel 


SCI 


Hiickel 


SCI 


0.00 


1.45 


3.26 


1.45 


3.26 


0.35 


1.44 


3.17 


1.42 


3.08 


0.70 


1.40 


3.13 


1.35 


3.01 


1.40 


1.28 


3.02 


1.14 


2.80 


1.80 


1.20 


2.94 


1.01 


2.65 



A. Lowest energy absorption. 

Our calculations of the optical transitions are within 
the SCI approximation. The SCI is known to give rea- 
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As seen in Table I, the optical gap decreases with t±. 
Interestingly, the decrease (relative to t± = 0) is greater 
for nonzero iJee than for H^e = 0, a result that is a con- 
sequence of increased delocalization of the excited elec- 
tron and hole of the 1B„ into the phenyl components 
for nonzero H^e, as we show below. The-|existing PPA 
and PDPA materials are finite oligomerstd, but for true 
long chain materials correct calculations of the optical ab- 
sorption should include self-consistent calculations of the 
bond alternation (instead of the rigid bond approxima- 
tion adopted here). In such a case, even larger decrease 
in the optical gap would be obtained. In the existing 
literature it is often assumed that PPA and PDPA have 
larger optical gaps than t-PA, but this is a consequence 
of the finite sizes of the substituted materials. The op- 
tical gaps of the substituted materials are smaller than 
the polyenes with the same chain length, with the gap in 
PDPA smaller than that in PPA. 

In Table II we have given the x- and y-components of 
the transition dipole couplings fj,x and iiy between the 
ground state and the 1B„ state of the PDPA oligomer 
(the electronic charge is taken as 1, and the dipole cou- 
plings are in units of A) . The transition dipole coupling is 
known to be a direct measure of the electron-hole corre- 
lation length in the IBu- We do not show similar results 
for PPA, as the difference is only quantitative. There 
are three interesting features of these results. First, for 
each t±, fix decreases for nonzero H^e, a result that is 
common to all other 7r-conjugated polymers, and is a 
signature of exciton formation and confinement of the 
electron-hole pairE^llZl in the 1B„ in the a;-direction for 
nonzero H^e- Second, t± decreases /i^ for both H^e — 
and ^ 0, but this decrease is considerably larger for 
nonzero Hee- For example, the decrease in fj,x between 
t± = and t± = 1.8 eV is 2.6 % for H^e = but 8.2 
% for Hee 7^ 0. This effect is reminescent of the ef- 
fect of t± on the optical gap (see above), and is a sig- 
nature of greater longitudinal electron-hole confinement 
for nonzero t±^ when Hee is nonzero. Note that with the 
Ohno Coulomb parameters, which are intermediate in 
strength, relatively small increase in the confinement can 
lead to 2Ag - 1B„ crossover. In the case of linear chain 
polyenes within the rigid bond approximation (i.e., hop- 
ping integrals to{l±S)) the difference. in the electron-hole 
correlation length in the 1B„ stateE^. between S = 0.07 
(E(2Ag) < E(1B„)) and S = 0.3 (E(2Ag) > E(1B„)) is 
only 7 %. In the present case, the increased longitudi- 
nal confinement is accompanied, as well as driven, by the 
electron correlation induced delocalization in the trans- 
verse direction, as observed from the very interesting be- 
havior of fly . As seen in Table II nonzero Hee iTLcreases 
fiy for arbitrary fixed which is completely contradic- 
tory to the behavior in the a;-direction. To the best of our 
knowledge, increasing transition dipole coupling with Hee 
has not been found in any other 7r-conjugated system. As 
noted in reference |l^ the fiy corresponding to an isolated 
trans-stilbene molecule oriented in the transverse direc- 
tion as in Fig. 1(b) decreases with Hee, thereby exhibiting 



the more common behavior. The opposite and surprising 
behavior in PDPA is a signature of greater electron-hole 
delocalization in the y-direction for the correlated IB^ 
than for the uncorrelated 1B„. Taken together, these 
results indicate a synergistic relationship between the 
confinement in the a;-direction and the enhanced delo- 
calization in the y-direction: the delocalization in the 
y-direction increases the confinement in the a;-direction 
and vice versa. 

The larger electron-hole delocalization in the y- 
direction is an electron correlation effect that can be 
understood from the natures of the one-electron levels. 
As shown in Fig. ||(a), the contribution of the back- 
bone chain atoms to the one-electron wavefunction of the 
highest occupied one-electron level (corresponding to the 
HOMO of the oligomer in the present context) decreases 
with t^. We have performed similar calculations for all 
one-electron levels, and for all the contribution of the 
chain atoms to the one-electron wavefunction decreases 
with increasing energy separation from the chemical po- 
tential. Now, for Hee = 0, the IB^ is simply the HOMO 
LUMO transition. For Hee 7^ 0, the 1B„ is a su- 
perposition of several transitions, including the HOMO 
-> LUMO, HOMO ~ 1 LUMO + 1, HOMO - 2 ^ 
LUMO -f 2, etc. The y-polarizations of these higher 
energy excited configurations are larger simply because 
the wavefunctions themselves have smaller contributions 
from the chain atoms and larger contributions from the 
phenyl group atoms, and as consequence the correlated 
1B„ wavefunction itself has a larger y-character. In Ref- 
erence ^5], we had given a configuration space argument 
for this electron correlation induced electron-hole delocal- 
ization in the y-direction. The present k-space argument 
is essentially the same explanation, except that the origin 
of the enhanced delocalization can be understood more 
quantitatively within these k-space arguments. 

As we have already pointed out in section III, increased 
confinement in the longitudinal direction can in princi- 
ple raise the energy of the 2Ag with respect to the IBu. 
Although this still needs direct verification, the impor- 
tant point here is that the increased confinement in the 
x-direction is a consequence of the increased transverse 
delocalization, which in turn reduces the optical gap (see 
Table I). 

TABLE II. The x and y components of lAg IBu tran- 
sition dipole moments fj.^ and fiy (in A), as a function of t±, 
for a PDPA oligomer with eight double bonds. Note t± — 0.0 
eV corresponds to an eight double-bond polyene. 



t± 




= 


Hee 


7^ (SCI) 




fix 




fJ.x 




0.00 


3.07 


0.35 


2.08 


0.73 


0.35 


3.06 


0.40 


1.97 


0.92 


0.70 


3.04 


0.54 


1.96 


0.99 


1.40 


3.00 


0.97 


1.93 


1.21 


1.80 


2.99 


1.24 


1.91 


1.34 
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FIG. 5. The calculated x (solid line) and y (dotted lines) components of the optical absorption spectra within the Coulomb 
correlated model for finite PPA (Figs, (a) and (c)) and PDPA (Figs, (b) and (d)) oligomers with eight double bonds on the 
backbone polyene chain. Note the difi'erent y-axis scales. In all cases t± = 1.4 eV, while to- = in Figs, (a) and (b) and 1.0 eV 
in Figs, (c) and (d). Coulomb interactions increase the j/-character of the lowest absorption and the a;-character of the higher 
energy absorption bands. 



B. Overall absorption spectrum 

In Figs, ^(a) and (b) we have shown the absorption 
spectra for PPA and PDPA for the case of nonzero Hcb 
but zero Hbb for the specific case of t± ^lA eV. We have 
included both x- and y-polarizations of the absorption 
bands in these figures. With nonzero H^e there occur 
absorption bands at energies as high as 10 eV. These 
are not included in the Figures. From comparison of 
these results with the absorption spectra of Figs. ^ (a) 
and (b), the lowest absorptions in both PPA and PDPA 
have much larger contributions from the y-component 
here than for Hee = 0, as has already been pointed out 
in the previous subsection. 

There is an even stronger difference between H^e — 
and Hee 7^ in the high energy region. For Hee = the 
dominant absorption in the high energy region is pre- 
dominantly y-polarized, although there is also substan- 
tial contribution from x-polarized transitions. In con- 
trast, for i/ge 7^ 0, the strong absorption bands centered 
at ~ 7 eV in both PPA and PDPA in Figs. 5(a) and 
(b) are not molecular absorptions but primarily charge- 
transfer chain benzene and benzene — > chain absorp- 
tions which have shifted to high energies because of e-e 
interaction. The one-electron chain-derived bands are 
predominantly delocalized along the a;-direction (with, 
however, substantial but smaller y-delocalization) and 
are also strongly hybridized with the benzene d, d* MOs. 
The high energy absorption bands here originate primar- 
ily from transitions between these delocalized levels and 



the benzene I, I* MOs, and as discussed in the context of 
one-electron theory, these transitions are x-polarized. 

The charge-transfer type absorptions are redshifted 
further upon inclusion of Hbb, as seen in Figs, ^(c) 
and (d) , where we have shown the optical absorptions for 
PPA and PDPA, respectively for =1.0 eV (note that 
the ^-contribution to the oscillator strength of the lowest 
energy absorption to the IB^ has remained unchanged, 
as was claimed in the previous subsection). The red- 
shift of the charge-transfer absorption also occurs within 
one-electron theory. However, within one-electron the- 
ory, there occurs a strong overlap between the a:-polarized 
charge-transfer absorption and the y-polarized molecular 
absorption (see Fig. |)(c) and (d)). In contrast, in the 
energy region of interest, there occurs only strongly a;- 
polarized charge-transfer absorptions, as the molecular 
absorptions are now strongly blueshifted by e-e interac- 
tions. The surprising effect of the e-e interactions is then 
that the absorption in the relatively high energy region is 
more strongly x-polarized than the lowest energy absorp- 
tion, even though the latter is predominantly a backbone 
chain absorption. 

Experimentally, a strong high energy absorption in 
PDPA derivative has already been seen by Tada et al. 
but this is in an unoriented sample and the polariza- 
tion characters are unknown. We predict that in ori- 
ented samples this absorption will be overwhelmingly x- 
polarized, which, as seen from comparisons of Figs. ^ 
and 1^, is strictly a correlation effect. In the absence of 
e-e interaction, the second band should have mixed x and 
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y-character. 



VI. 2A„—lBu ORDERING 



In our earlier Letter, we presented preliminary cal- 
culations on the relative ordering of 2Ag and ex- 
cited states in PDPAs, using selected subsets of MOs 
of a two jdDuble-bond oligomer, within the fuU-CI (FCI) 
approachEj. Here we present more sophisticated calcu- 
lations using all the MOs of oligomers of PDPAs and 
PPAs with four as well as six backbone carbon atoms, 
by employing the MRSDCI approaclJfl. The MRSDCI is 
an efficient size-consistent algorithm that has been shown 
to give energies of the 1B„ and the 2Ag that are as ac- 
curate as the quadruple-CI (QCI) approach for polyenes 
up to sixteen carbon atomsQ, and has also been used 
to determine the frequency depjsndent nonlinear opti- 
cal properties of short polyenea^j. Note that the PDPA 
oligomer with six backbone carbon atoms contains 42 
carbon atoms altogether, and the Hamiltonian matrices 
are much larger than those previously investigated. As 
shown below, strong emphasis was placed on obtaining 
very high accuracy, and the oligomer with six backbone 
carbons is the largest that could be investigated without 
loss of accuracy. 

Our calculations are within the rigid bond approxima- 
tion, i.e., the hopping integrals along the chain are taken 
to be to (1 ±(5). We know that enhanced bond alternation 
can lead to 2Ag state being higher than the IBu statdj'Q. 
Therefore, in order to understand the influences of elec- 
tron correlations and bond alternation on the 2Ag-\Bu 
ordering separately, we considered backbone geometries 
of PPAs and PDPAs both with standard polyene bond 
alternation [5 = 0.07), and without any bond alternation 
(i5 = 0.0). For the case of PDPA oligomers, full use of the 
Ci symmetry of the system was made during the corre- 
lated calculations. Therefore, Ag andJ3„ subspaces were 
diagonalized in separate calculationsO. However, for the 
PPA oligomers, because of the lack of any symmetry, Ag 
and Bu states were obtained by diagonalizing the same 
Hamiltonian matrix. The \Bu state was identified in this 
case as the first excited state with large dipole coupling 
with the ground state. In the MRSDCI calculations pre- 
sented below, we have ignored H^b, based on our results 
in Figs. 3, 4 and 5 which clearly indicate that Hbb plays 
an insignificant role in the low energy physics. As far as 
the phenyl-chain interaction term Hcb is concerned, in 
all the calculations reported below it was parametrized 
with t^^lA eV. 

The methodology behind the MRSDCI calculations 
is as follows. The calculations are initiated with a 
Hartree-Fock (HF) computation of the ground state of 
the oligomer concerned, followed by a transformation of 
the Hamiltonian from the site representation to the HF 
MO representation. Subsequently, a singles-doubles CI 
(SDCI) calculation is performed, the three states of in- 



terest, viz., the \Ag, the 2Ag^ and the 1B„ are examined, 
and the N^e/ configuration state functions (CSFs) mak- 
ing significant contributions to their many-particle wave 
functions are identified. The next step is the MRSDCI 
calculation for which the reference space consists of the 
Nre/ CSFs identified in the previous step, and the over- 
all Hamiltonian matrix now includes configurations that 
are doubly excited with respect to these reference CSFs 
(thereby including the dominant quadruply excited con- 
figurations). The new lAg, 2Ag, and IBu states are now 
re-examined to identify new CSFs contributing signifi- 
cantly to them so as to augment the reference space for 
the next set of MRSDCI calculations. This procedure is 
repeated until satisfactory convergence in the excitation 
energies of the 2Ag and IB^ energies is achieved. By 
the time convergence is achieved typically all CSFs with 
coefficients of magnitude 0.01 or more in the correspond- 
ing many-particle wave functions have been included in 
the MRSDCI reference space. Naturally, this leads to 
very large CI matrices, e.g., the dimension of the CI ma- 
trix in the MRSDCI calculation on the B^ space of the 
three unit PDPA oligomer with Nref = 12, was « 2.22 
millions. The lowest few eigenvalue and eigenvectors of 
such la*ee matrices were obtained by the Davidson pro- 
ccdurco. The very large N^e/ required for convergence 
(even for the IB^,) is direct evidence of the contributions 
by the benzene MOs to the correlated wavefunctions of 
PPA and PDPA. 

The convergence pattern of our MRSDCI calculations 
with respect to N^e/ is demonstrated for the case of the 
PPA oligomer with four backbone carbon atoms in Ta- 
3. These calculations correspond to (5 = 0. The 
g and IBu excitation energies of various intermedi- 
ate MRSDCI calculations are measured with respect to 
the best lAg ground state obtained in the largest MRS- 
DCI calculation the table {Nref — 25). It is for this 
reason that with increasing Nref, the 2Ag and IB^ ex- 
citation energies decrease monotonically. It is clear from 
the behavior of excitation energies with respect to Nref 
depicted in Table |l| that by the time the MRSDCI cal- 
culation with the largest value of N^e/ is performed, the 
2Ag and 1B„ excitation energies have converged to an 
acceptable accuracy. 

Our overall MRSDCI results for 2Ag and 1B„ ener- 
gies for four and six backbone atom PPA and PDPA 
oligomers (both with and without bond alternation), 
along with the corresponding final N^e/'s, are summa- 
rized in Table IV. In the same Table, for the sake of 
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comparison, the results of full CI calculations on four 
and six atom polyenes, both with and without bond al- 
ternation, are also pre sented. Upon examining the data 
presented in Table the following trends emerge: (i) 
in PDPA oligomers without backbone bond alternation, 
IBu and 2Ag states are nearly degenerate, while, in the 
same oligomers with bond alternation, the IBu state is 
significantly below the 2Ag, (ii) in PPA oligomers with- 
out bond alternation the 2Ag state is below the IB^, 
while for oligomer geometries with bond alternation 2Ag 
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is very slightly above the 1B„, (iii) in all cases, the dom- 
inant effect of the phenyl-substitution is the lowering of 
E(1B„), which decreases as t-PA > PPA > PDPA. This 
last result is in agreement with the SCI results for much 
longer oligomers in Table ^ confirming our statement 
that the SCI captures the essential nature of the 1B„ 
state. 

The occurrence of the 2Ag slightly below the IBu in 
the PDPA oligomer with S — and six backbone carbons 
is not surprising. As discussed in reference exactly as 
bond-substituted systems can be thought of as effective 
linear chains with large bond alternation, site-substituted 
materials can be thought of also as linear chains, but now 
with composite sites (as opposed to bonds) with small ef- 
fective on-site Coulomb repulsion Ugff- In linear chains 
of sufficient length, E(2Ag) < E(1B„) for 5 = for any 
Ueff- The more important point, however, is that the 
energy gap E(1B„) - E(2Ag) has decreased from 1.15 eV 
to 0.02 eV in going from the polyene with S = to the 
PDPA oligomer. As discussed in section IV. A, the per- 
sistence of the Fermi surface degeneracy for 6 = 0, even 
upon phenyl-substitution (see Figs. 2(a) and (b)), ensures 
that unconditional bond alternation occurs in both PPA 
and PDPA, and therefore the nonzero 6 results in Table 
IV are more appropriate for the real systems. Although 
our calculations are for relatively short oligomers, we ex- 
pect our conclusions regarding the 2Ag — IBu ordering 
in PDPAs to be valid even for long chains (in particu- 
lar, for the experimental systems, which are also finite, 
albeit of length that is twice the lengths studied here). 
To begin with, the size-dependence of the 2Ag - 1B„ en- 
ergy gap, for the case E(2Ag) > E(1B„) is steepest at 
the smallest sizes (a result easily confirmed within one- 
electron theory). Nevertheless, E(2Ag) - E(lBtj) has al- 
most converged already for the PDPA oligomers with 6 ^ 
0: for the oligomer with four backbone carbon atoms this 
energy difference is 0.4 eV, while for the oligomer with 
six backbone carbons the difference is 0.37 eV. Second, 
our calculations are for rather large Coulomb interactions 
and fairly small t± , while in the ceal systems the Coulomb 
interactions are perhaps smallerEZI and t± might be larger. 
Finally, it is clear from Table tV that even if the 2Ag 
were to be below the 1B„ in short PDPA oligomers, the 
energy difference would be tiny. Density matrix renor- 
malization group calculations for linear chain systems 
have suggested that specifically in these cases there is 
a size-dependent 2Ag - 1B„ nrasspver, where in longer 
chains E(2Ag) > E(1B„) againEaO. Thus from aU possi- 
ble considerations, the lowest two-photon state is above 
the lowest one-photon state in PDPA. In contrast to the 
PDPA oHgomers, E(2Ag) - E(1B„) has not converged 
in the PPA oligomers, where this number is 0.13 eV for 
the oligomer with four backbone carbons and 0.06 eV in 
the oligomer with six backbone carbons. It is conceiv- 
able therefore that in slightly longer oligomers, the 2Ag 
is slightly below the IB^ in PPA. 



VII. CONCLUDING DISCUSSIONS 

We begin this section with the summary of results ob- 
tained in this paper. We agree with the authors of the 
experimental reference |l^ that light emission in PDPAs 
involves predominantly the backbone chain and not the 
trans-stilbene moiety. However, as we have shown here, 
the lowest absorption has a strong y-character, which is 
enhanced by e-e interactions, and it is this enhanced delo- 
calization in the y-direction that is ultimately the origin 
of the 2Ag occurring above the IBu in these systems. 

TABLE III. Convergence pattern of MRSDCI 2Ag and 
IBu energies for PPA with four backbone carbon atoms and 
uniform hopping integrals with respect to increasing N^ef- 
Nref denotes the total number of reference configurations in- 
cluded in the MRSDCI calculation in question. The 2Ag and 
IBu energies are with respect to the lAg ground state ob- 
tained in the largest calculation {Nref = 25). AU energies are 
in eVs. 

System 
^PA 



Nref 


E{2Ag) 


E{lBu) 


2 


5.50 


4.77 


6 


4.94 


4.69 


8 


4.74 


4.65 


10 


4.67 


4.65 


12 


4.60 


4.63 


15 


4.52 


4.62 


19 


4.50 


4.61 


21 


4.46 


4.60 


23 


4.43 


4.60 


25 


4.42 


4.59 



TABLE IV. The best MRSDCI results for 2Ag and IBu 
energies of the oligomers of PDPA and PPA containing four 
and six carbon atoms along the polyene backbone. Backbone 
configurations with both S = 0.0 and S — 0.07, were consid- 
ered. For the sake of comparison, the results of the full CI 
calculations (indicated by FCI in the A'^,.e/ column) for the 
two states of the corresponding polyenes are also given. Since 
for PPA oligomers the Ag and Bu states were obtained in the 
same calculations, the same values of Nref are indicated for 



both the Ag and Bu states. 



Energy Gaps (eV) 
Four backbone carbons Six backbone carbons 
System State Nref 5 = 0.0 5 = 0.07 Nref S = 0.0 S = 0.07 



PDPA 


2Ag 


17 


4.38 


4.98 


13 


3.62 


4.39 




IBu 


8 


4.25 


4.58 


12 


3.64 


4.02 


PPA 


2Ag 


25 


4.42 


5.03 


14 


3.66 


4.47 




IBu 


25 


4.59 


4.90 


14 


4.02 


4.41 


polyene 


2Ag 


FCI 


4.66 


5.35 


FCI 


3.51 


4.37 




IBu 


FCI 


5.49 


5.80 


FCI 


4.66 


5.02 
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The authors of reference emphasize exciton confine- 
ment due to relaxations involving phonons. As shown 
here, enhanced exciton confinement along the backbone 
is a consequence of the electron-hole dclocalization in the 
transverse direction. This does not imply that electron- 
phonon interactions are unimportant. Rather, even if 
relaxations involving phonons play a strong role in the 
experimental systems, the 2Ag had to be already close 
to the lHu in the PDPAs (given the large separation be- 
tween these in t-PA) due to a purely electronic mecha- 
nism, and the present work gives the details of this elec- 
tronic mechanism. We emphasize that the longitudinal 
confinement does not necessarily imply enhanced exciton 
binding energy, as the dclocalization in the transverse 
direction will also reduce the energy of the conduction 
band threshold. The bond alternation along the back- 
bone decreases with phenyl substitution, also as a direct 
consequence of smaller chain atom contributions to the 
one-electron levels at and near the Fermi level. Inclusion 
of the overlaps between the benzene MOs that occur due 
to phenyl ring torsions does not change the results of any 
calculation that probes the low energy photophysics. Fi- 
nally, we predict a strong second high energy absorption 
that will be predominantly x-polarized in oriented ma- 
terials, in contradiction to the prediction of one-electron 
theory. Experimental work can therefore directly verify 
the role of e-e interactions. 

In the case of PPAs, our results are more difficult to in- 
terpret. The convergence behavior of E(2Ag) andE(lB„) 
in the two oligomers that have been studied suggests that 
here E(2Ag) < E(1B„) in longer chains, although the en- 
ergy difference can be very small. The small energy dif- 
ference between these levels would imply that even small 
perturbations can either reverse this ordering, or lead 
to strong state mixing. We speculate that the behavior 
of the PAPAs is related to this proximity between the 
2Ag and the 1B„. PAPAs with small alkyl groups prob- 
ably are similar to the PPAs of the same length and are 
therefore weak emitters. Large alkyl groups either lead to 
greater main chain bending (and therefore greater exci- 
ton confinement) and state reordering, or because of their 
larger electron donating capability they cause a greater 
electronic perturbation leading to greater state mixing 
between the 2Ag and the 1B„. 

In our investigation of optical properties here, we have 
focused on ground state absorption and theoretical 
studies of excited state absorptionstllEa is a topic of fu- 
ture research. In the case of PPP and PPV, it has been 
claimed that the occurrence of multiple classes of one- 
electron bands leads to several different photoinduced ab- 
sorptions in ultrafast pump-probe spectroscopyEj, and it 
is conceivable that similar results will be found here too. 

We make one final remark regarding the future ap- 
plication of this work for the design of novel materials. 
Photo- and electroluminescence of 7r-conjugated poly- 
mers are currently the focus of intense research, and re- 
search on the design of novel light emitting diodes usiiic 
these materials has already reached a mature stagec^l. 



Initial work has also shown proraiaina; results for laser ap- 
plications using these materialsETcj. In all these cases, 
however, the materials being used correspond to those 
obtained by bond substitution, and the optical gaps of 
such systems are necessarily larger than that of t-PA. Our 
demonstration that light emission can be a consequence 
of enhanced dclocalization in a direction transverse to the 
direction of extended conjugation, and that the optical 
gaps in this latter class of systems are smaller than linear 
polyenes of the same length, introduces the interesting 
possibility of designing novel rigid long chain light emit- 
ting TT-conjugated polymers with optical gaps smaller 
than that in t-PA. The actual synthesis of rigid mate- 
rials may be difficult, as the finite conjugation lengths of 
PPAs and PDPAs already demonstrate, but nevertheless, 
the present work shows this to be conceptually feasible. 
Infrared dye lasers are extremely rare, and the possibil- 
ity of obtaining new light sources at telecommunications 
wavelengths and in the eye-safe spectral region should 
make this an exciting new research direction. 
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